Probing the Leggett-Garg Inequality for Oscillating Neutral Kaons and Neutrinos 



o 

(N 

u 

Oh 
< 
OS 



43 

Oh 
i 

cr 



(N 

O 
m 

> 

'i—i 

X 



D. Gangopadhyay 1 (debashis@rkmvu.ac.in), D. Home 2 

(dhome@bosemain.boseinst.ac.in) and A. Sinha Roy 1 (animesh@rkmvu.ac.in) 

1 Department of Physics, Ramakrishna Mission Vivekananda University, Belur Matth, Howrah-7 11202, India and 

2 CAPSS, Department of Physics, Bose Institute, Sector-V, Salt Lake, Kolkata 700 091, India 

The Leggett-Garg inequality (LGI) as a temporal analogue of Bell's inequality, derived using the 
notion of realism, is applied in a hitherto unexplored context involving the weak interaction induced 
two-state oscillations of decaying neutral kaons and neutrinos. The maximum violation of LGI 
obtained from the quantum mechanical results is significantly higher for the oscillating neutrinos 
compared to that for the kaons. Interestingly, the effect of CP non-invariance for the kaon oscillation 
is to enhance this violation while, for neutrinos, it is sensitive to the value of the mixing parameter. 
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Introduction-Nonclassical features of the quantum 
world and their deeper implications have been explored 
since the inception of quantum mechanics (QM). A novel 
slant on this line of study was provided by Bell's inequal- 
ity (BI)[T] whose distinctiveness stems from the feature 
that it is a testable algebraic consequence of a very basic 
notion known as local realism. This assumes that ob- 
scrvables pertaining to any object, even when not mea- 
sured, have definite values (the notion of realism), and 
that results of individual measurements of these prop- 
erties remain unaffected by spatially distant events (the 
locality condition). BI essentially sets a bound on a cer- 
tain combination of correlation functions corresponding 
to outcomes of measurements on two spatially separated 
systems. For suitable relative orientations of these mea- 
surements, BI is violated by the relevant QM results 
for appropriate states of the entangled systems. Exten- 
sive experimental investigations [2] over the past three 
decades have succeeded in closing all possible loopholes, 
thereby providing convincing empirical repudiation of BI, 
consistent with the QM predictions. The upshot of these 
studies is the realization that any realist description of 
microphysicone of the earliest examples of CP violation 
in natureal phenomenon has to be intrinsically nonlocal. 

A stimulating twist to the above line of study was pro- 
vided by the Leggett-Garg inequality (LGI) [3J 2] which 
is a temporal analogue of BI in terms of time-separated 
correlation functions corresponding to successive mea- 
surement outcomes for a system whose state evolves in 
time. Notion of realism is invoked in deriving LGI by 
assuming that a system, during its time evolution, is at 
any given time in a definite one of the available states. 
Instead of the locality condition, the notion of noninva- 
sive measurability (NIM) is used. This means assuming 
that it is possible, in principle, to determine which of the 
states the system is in, without affecting the state itself 
or the system's subsequent dynamics. Leggett [3] has 
argued justifying why NIM is to be considered a "nat- 
ural corollary" of the notion of realism. QM violation 
of LGI in suitable examples would, therefore, signify re- 
pudiation of the notion of realism that includes the as- 



sumption of NIM. Thus, while furnishing a signature of 
distinctly quantum behaviour, LGI can be regarded as 
complementing BI in providing valuable insight into the 
nature of physical reality that is entailed by nonclassi- 
cality of quantum systems. Hence it has been of consid- 
erable interest to investigate the extent to which LGI is 
violated by QM for various types of systems. The orig- 
inal motivation [3j |4] that led to LGI was to use it for 
probing the limits of quantum mechanics in the macro- 
scopic regime; e.g., in the context of suitable experiments 
involving the rf- SQUID device [5]- In recent years, study 
of QM violation of LGI has been carried out for a vari- 
ety of solid-state qubit systems [6] and optical systems 
7\. Empirical violations of LGI have been shown using 
'weak measurements' [5], employing liquid-state nuclear 
magnetic resonance in chloroform |9 , as well as for sin- 
gle spins in a diamond defect center |10j and for nuclear 
spins processing in an external magnetic field [I II . 

Against this backdrop, the present paper initiates stud- 
ies along an earlier unexplored direction. Neutral kaons 
and neutrinos are particularly interesting systems in el- 
ementary particle physics which provide remarkable ex- 
amples of oscillations in time over different states - an 
inherent property of these systems which is governed 
by the fundamental electro-weak interaction. While the 
phenomenon of neutrino oscillation |12j has played a vi- 
tal role in establishing non-zero rest mass of neutrinos, 
the significance of kaon oscillation [13] is that it involves 
eigenstates of the effective weak interaction Hamiltonian 
which exhibit decay embodying CP non-invariance - the 
earliest detected example of CP violation in nature. For 
these systems, it is, therefore, of prime interest to investi- 
gate not only the extent to which QM results can violate 
LGI, but also the nature of dependence of QM violation 
of LGI on certain relevant key parameters, such as the 
magnitude of CP violation for kaon oscillation and the 
mixing angle for neutrino oscillation. Here it is relevant 
to note that an incompatibility between BI and the rele- 
vant QM results for entangled pairs of neutral kaons has 
been demonstrated using Wigner's version of Bell's in- 
equality [2], while there were also earlier attempts [TS] 



towards demonstrating QM violation of local realism us- 
ing neutral kaons. In view of such studies, it is thus of 
added relevance to examine the QM violation of a tem- 
poral analogue of BI for oscillating neutral kaons. 

The Leggett-Garg inequality - We begin by setting up 
the relevant form of LGI that will be used in this paper. 
For this, we focus on a two-state system whose temporal 
evolution consists of oscillations between the states, say, 
1 and 2. Let Q(t) be an observable quantity such that, 
whenever measured, it is found to take a value +1(— 1) 
depending on whether the system is in the state 1(2). 
Next, consider a collection of runs starting from identi- 
cal initial conditions such that on the first series of runs 
Q is measured at times t\ and £2, on the second at £2 and 
£3, on the third at £3 and £4, and on the fourth at £1 and 
£4 (here £1 < £2 < £3 < £4). From such measurements, 
it is straightforward to determine the temporal correla- 
tions Cij = (Q(ti)Q(tj)). Now, as argued by Leggett 
and Garg [31 H], it is possible to adapt in this context, 
the standard argument leading to a Bell-type inequal- 
ity with the times ti playing the role of apparatus set- 
tings. One can then use the following consequence of the 
assumptions of realism and NIM that were mentioned 
earlier. For any set of runs corresponding to the same 
initial state, any individual Q(U) has the same definite 
value, irrespective of the pair Q(U)Q(tj) in which it oc- 
curs; i.e., the value of Q{t{) in any pair does not depend 
on whether any prior or subsequent measurement has 
been made on the system. Consequently, the combiaar 
tion[Q{ti)Q{t2) + Q{t 2 )Q{h) + Q{h)Q{U)-Q{h)Q{U)] 
is always +2 or —2. If all the individual product terms 
in this expression are replaced by their averages over the 
entire ensemble of such sets of runs, the following form 
of LGI is then obtained 
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The above is, thus, an inequality imposing realist con- 
straints on the time-separated joint probabilities pertain- 
ing to oscillations in any two-state system. We now anal- 
yse its incompatibility with QM in the specific cases of 
kaon and neutrino oscillations respectively. 

Oscillating neutral kaons - Neutral kaons K a ,K° are 
pseudo scalar mesons, each of which is the antiparticle of 
the other and is distinguished by the strangeness quan- 
tum number S = +1 (-1) for K° (K°). Evolving un- 
der the effective weak interaction CP violating Hamilto- 
nian H — M — iT/2 where M and T are the Hermitian 
mass and decay matrices respectively, K° and K° de- 
cay, while giving rise to K°,K° oscillations. Eigenstates 
of H are respectively the long-lived (rj, ~ 10~ 7 s) and 
short-lived (rg ~ 10 _10 s)states \Kl,s) with eigenvalues 
^l,s = to l,s — *7l,s/2 corresponding to masses i~iil,s 
and characteristic decay widths 7l,s(= ^/tl.s) where 
\KL,s(t)) — exp(—i\L t st)\KL_s)- In terms of mutually 



orthogonal strangeness eigenstates K°,K° one has 
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where e is the complex CP violating parameter and 
(K L \K S ) = 2i?e(e)/(l + |e| 2 ). Using Eq. ([2} , and the time 
evolution of \Kl^s), one can calculate the time evolution 
of an initial pure K° beam. The probabilities of finding 
K° and K° states after time t are respectively 

., e" 7Lt + e- 7s * + 2e- 7 *cos(Am£) ,, 

p K° (t) = -; (3) 



^°W-|(i + e )| 2 



e -7i.< 4. e -7s* _ 2e" 7 * cos(Am£) 
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where Am = m^ — TO5 and 7 — (75 + 7^)/2. The joint 
probability P K a ^o(£i,£2) of finding the states \K°), \K°) 
at the respective times £1 and £2 is then 
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Similarly, one can calculate the other joint probabilities 
PK°,K°(ti,t2),PK0^ K a(t 1 ,t 2 ) and Pjco i Ko(ti,t2). Then 
the time correlation function Ci2(= \Q(ti)Q(t2)}) can 
be evaluated by using expressions like Eq. ^ to obtain 
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+Re(e)e~'> tl {e-'< L< - t2 - tl ) + e -7s(*2-ti) j. cos (Am£i) 
(i+| e | 2 )( e -^'i+ e -"'s t i) r c ^^ f ,( t ^t 1 ) | e - 7s ( t3 _ tl )j 



+2i?e(e)e- 7 * 2 cos(Am£i) cos{Am(£ 2 - £1)} (6) 

An interesting point to be noted here is that, in the pres- 
ence of CP violation (e 7^ 0), C12 displays dependence 
on both the temporal separation (£2 — £1) and on t\. If 
e = 0, the dependence on t\ disappears, as can be checked 
from Eq.([6]). The other temporal correlation functions 
C23, C34 and C14 can be calculated in the same way and 
they also show a similar feature. Here we may note that 
under CPT invariance, CP violation implies time reversal 
non-invariance. This explains why in the presence of CP 
violation, the C'^s are not merely functions of the tem- 
poral separation but also depend on the initial instant of 
the pertinent time interval. 

Next, using equations like Eq.([6]), one can evaluate C 
as defined in Eq. ([I]) in order to study the QM violation of 
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FIG. 1: In Fig. 1(a) pertaining to kaon oscillation in the pres- 
ence of CP violation, the overall behaviour of the QM cal- 
culated value of C as a function of (At/rs) is shown within 
the oscillation time range, where we take ti/rs — 5.3. In 
Fig. 1(b), we focus around the region where the QM value of 
C is maximum as (At/rs) is varied. Comparing the unbroken 
and broken curves for the cases with and without CP viola- 
tion, one can see the enhancement of the QM violation of 
LCI due to CP violation. In Fig. 1(c), pertaining to neutrino 
oscillation, we show the overall behaviour of the QM value 
of C as the parameter AL/E is varied for the KamLAND 
experimental setup. 



LGI for kaon oscillations. By varying the choices of the 
time intervals, it is found that the maximum value of C is 
attained essentially when the temporal intervals are cho- 
sen to be the same, i.e., when £4 — ^3 = £3 — £2 = £2 — ^i = 
At. The detailed expression for C under this condition 
and in the presence of CP violation is given in the sup- 
plementary material[l6 where it can be seen explicitly 
that C depends on both At and t\. The experimen- 
tally determined values of Am, r s , tl, |e| and Re(e) occur- 
ring in the expression for C are r s = 0.8958 x 10 _10 sec, 
r L = 0.5084 x 10~ 7 sec, Re(e) = 1.596 x 10~ 3 [T7], while 
Am = 3.843 x 10~ 12 MeV, |e| = 2.232 x 10~ 3 [IS]. We 
consider various choices of t\ and At within an appropri- 
ate time scale (ti/rs < 10, At/rs ^ 10) over which the 
decaying kaons show appreciable oscillation. It is then 
found that the QM calculated maximum value of C in the 
oscillation time range is 2.36463 when At = 0.7889 x rs 
and ti — 5.3 x 7-5. Note that the maximum QM viola- 
tion of LGI in this case is significantly smaller than the 
upper bound [19] of the QM violation of the Bell-type 
inequality given by 2\/2 = 2.82843. For a given value of 
ti = 5.3 X rg, the variation of QM calculated quantity C 
with At/rs is shown as the unbroken curve in Fig. 1(a). 
Here we observe that, if, for example, the difference be- 
tween the decay parameters (jl — Tg) had a lower value, 
the maximum QM value of C would have increased. 

If we ignore the effect of CP violation and put e = 0, 
C depends only on At as can be seen from the detailed 
expression for C given in the supplementary material |16j. 



By varying At, it is then found that the maximum QM 
value of C in the oscillation time range is 2.36448 for 
At = 0.7890 x rs, whereas, as mentioned earlier, the 
maximum value of C in the presence of CP violation is 
2.36463. Thus, while the overall behaviour of the QM cal- 
culated quantity C as At/rs is varied is the same with 
and without CP violation, a notable feature is that the 
QM violation of LGI is enhanced by an amount 0.00015 
in the presence of CP violation (see Fig. 1(b)). We 
have also checked that if |e| was, for example, larger, say, 
2.23 x 10 -2 , the maximum QM value of C would have 
increased to 2.36667. Conceptually, this means that the 
non-invariance of CP symmetry in weak interaction re- 
sults in an enhanced non-classicality of the kaon oscilla- 
tion as shown by an increased QM violation of LGI. 

Oscillating neutrinos - It has been well established 
by a number of experiments that during propagation, 
neutrinos undergo oscillations among three flavor eigen- 
states i/ e , z-^ , ^ r . Here our treatment of LGI using the 
neutrino oscillation phenomenon is analysed in the con- 
text of the KamLAND experimental setup where the 
effect of the mixing angle #13 can be considered neg- 
ligibly small |20j . Then for neutrinos, one can essen- 
tially consider two-flavor oscillation involving transitions 
between v e and v^. Here the mass eigenstates are Vi 
and V2 with energy eigenvalues E12 respectively, where 
Vi = cos Qv^ — sin 9v e , vi = sin 91/^+ cos 9v e and ^1,2 (*) = 
Vi^exp^—iEi^t/h). Now, given an initial electron neu- 
trino beam, after time t, the probability of obtaining v^ 
and v e are respectively [21] 
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where 9 is the mixing angle, Am = m\ — ra\ where 
mi and m,2 are masses corresponding to the states Vi 
and i>2 respectively and E is the mean energy of neu- 
trino mass eigenstates. Next, one can calculate the four 
joint probabilities (similar to those defined for the kaon 
oscillation) P Ve ,^(ti,t2)»-Pi/„i/«(*i,*2),iV (1 ,i/ e (*i)*2) and 
Pf ,u (^l)^)- The time correlation function C\% in this 
case is given by 



C12 — 
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c 4 (t 2 -h) 



(9) 



Similarly, the functions C23,C34 and C14 can be calcu- 
lated. Here the temporal correlation functions depend 
only on the temporal separation. Next, using equations 
like Eq.Q, one can evaluate the quantity C as defined in 
Eq. |l]) in order to study the QM violation of LGI for the 
oscillating neutrinos. By varying the choices of the time 
intervals, similar to the case of kaon oscillation, it is found 



that the maximum QM value of C for the neutrino oscilla- 
tion, too, occurs only when £4—^3 = £3— ti — £2—^1 = AL 
The quantity C under this condition becomes a function 
of At/E. Now, if AL be the distance traversed by the 
neutrinos during time interval At, then C is given by 



C = 
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(10) 



with AL — cAt. Since neutrinos are non-decaying, it 



is clear from Eq.(10) that, unlike in the case for neutral 



kaons, the QM calculated quantity C continues to oscil- 
late with time with maxima points occurring at various 
points of (AL/E) . The experimentally obtained values 
of Am 2 c 4 , tan 2 8 in the case of the KamLAND experi- 
mental setup (involving reactor neutrinos and geo neu- 
trinos) are respectively 7.58 x 10 -5 ev 2 and 0.56 [5D]. For 
such data, the maximum QM value of C is 2.76 and this 
maximum value is repeated at various points of (AL/E); 
see Fig. 1(c). Further, we have checked by varying Am 2 
over a range of values that there is no significant change 
from the maximum QM value of C. On the other hand, 
the maximum QM value of C depends sensitively on the 
value of the mixing angle 8. If 8 had a higher value 
than the experimentally determined value, the maximum 
QM value of C would have increased, reaching the upper 
bound 2\/2 for 8 = tt/4, and then would have decreased 
if the value of 8 was still higher. 

Concluding remarks - Comparing the magnitudes of 
the QM violation of LGI in the cases of the kaon and two- 
flavour neutrino oscillation, it is found that while in the 
former case, the ratio by which the realist upper bound 
given by LGI is violated by the relevant QM results is 
18%, this ratio is 38% for the neutrino oscillation, thereby 
significantly higher than the corresponding value for the 
kaon oscillation, but smaller than the upper bound of this 
ratio (i.e., 41%) by which LGI can be violated by QM. 
On the other hand, for the kaon oscillation, the effect of 
CP violation is to enhance the QM violation of LGI by 
0.008% - albeit small, but not a negligible effect that is 
conceptually interesting. 

Possible implications of the above results call for care- 
ful reflection. Besides, it could be worthwhile to make a 
detailed quantitative comparison of the results obtained 
in this paper with the QM incompatibility with local re- 
alism analysed for neutral kaons O S3- Finally, we 
should mention that the phenomenon of neutrino oscil- 
lation recently studied in the context of the Daya Bay 
reactor neutrino experimental setup ; 22J involves a non- 
negligible value of the mixing angle #13, thereby entail- 
ing three-flavour neutrino oscillation. It could, therefore, 
be interesting to examine whether the QM incompatibil- 
ity with the notion of realism underpinning LGI can be 



probed in such a context as well. 
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